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HEAT TRANSFER TO CLOUDS OF FALLING PARTICLES 


I. InrRopucTION 


1. Object of Work.—This investigation was undertaken for the 
purpose of obtaining information essential to the design of equipment 
for heating finely divided materials falling through a gas. The im- 
mediate object of the work was the design of a flash calciner for the 
decomposition of hydrated zine sulphite, the product of one operation 
in a cyclic process for the recovery of sulphur dioxide from waste 
gases. A summary of the results and their application to that par- 
ticular problem have been given in a previous bulletin.* Because of 
the general applicability of the work, the results are presented here in 
more detail along with the theoretical treatment essential to a complete 
development of the problem. 

There is a growing interest in space reactions of solids, or between 
solids and gases, at high temperatures, because of the high rates of 
heat transfer to clouds of solid particles suspended in a gas. The 
method was suggested at one time for the carbonization of pulverized 
coal to produce “bubble coke.”+ Adherence of the particles to the 
retort walls, however, was not entirely prevented and the method did 
not gain favor, although it was recognized that the rate of heat trans- 
fer was rapid, and that conduction through the walls was actually a 
limiting factor on the capacity.t Application to the briquetting of 
blast furnace dust, the smelting of copper and zinc ores and for the 
production of soluble phosphates has been suggested.§ Other possible 
applications include the decomposition of sodium bicarbonate and the 
dehydration of clays and minerals. Several of these have been in- 
vestigated in the experimental equipment to be described. 

The general problem is also of interest in the design of spray 
dryers. Lapple and Shepherd{ have recently made a theoretical study 
of the trajectories of the particles in such equipment. From the infor- 
mation given in their paper, it should be possible to apply the basic 
equations developed here for heat transfer to determine the heat re- 
moved from the particles and, thus, the point at which solidification, 
or crystallization is complete. 

Anse CHSC Beets Wie ian. Cont. 08 Bit. Coal, Carnegie Inst., p. 419, 1926; W. 
Runge, Ibid., p. 697; F. M. Gentry, Combustion, vol. 20, p. 225, 1929. 

tF. B. Hobart and D. J. Demorest, Ohio State Univ. Eng. Exp. Sta. Bul. 65, 1932. 


SE. GC. St. Jacques, Ind. Eng. Chem. News Ed., vol. 15, p. 29, 1937. 4 
{C. E. Lapple and C. B. Shepherd, Ind. Eng. Chem., vol. 32, p. 605, 1940. 
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Il. THEroreticaL Discussion 


3. Nomenclature —Throughout this bulletin the equations are 
written in a form in which any consistent system of units may be 
used. In those cases in which designation of the units has been neces- 
sary, the English system has been used for quantities having immedi- 
ate practical significance, and the metric system for quantities used 
principally in the laboratory. This policy, rather than one of strict 
uniformity, was adopted in order to secure greater clarity for the 
average reader. 


a = average projected area of a single particle 
A = surface area; subscript o refers to the total area emitting 
radiation, p to that of each particle, w to that of the in- 
ternal wall of the furnace 
b, B = arbitrary constants 
c = specific heat of gas 
C = concentration of particles at any point in the cloud, i.e 
the number per unit volume 
Cp = drag coefficient, defined by Equation (60) 
d = differential operator 
D = internal diameter of furnace tube 
D, = arithmetic average diameter of particles in the cloud 
i 
F 


e) 


= fraction of total area A covered by dises 
= feed rate, in general, expressed as lb./(min.) (sq. ft.) 
rn = angle-emissivity factor for net radiation transfer 
g = acceleration due to gravity 
h = heat transfer coefficient; subscripts have the following 
meanings: m overall coefficient, r radiation coefficient, 
cp convection coefficient for spherical particle surface, ew 
convection coefficient at furnace wall. All coefficients ex- 
cept hw are based on particle area 


if = 


J, = 
k= 
= 
l= 


Aot = 


T = 


U,v, WwW = 


Qi = 
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intensity of radiation, subscript 1 refers to the intensity at 
point / measured from the front of the cloud, o refers to 
that at which / = 0, n refers to that in the direction normal 
to the surface 

Bessel function of the first order and the first kind 
thermal conductivity 

volume-shape factor, defined by Equation (59) 

length of a radiant beam; subscript o refers to the total 
length of beam 

length of path of particle, i.e., in general, the vertical dis- 
tance along the furnace 

a constant, defined by Equation (16) 

number of particles present in the furnace at any instance 
rate of heat transfer ; 

radius; subscript / refers to radius of a spherical particle, 
subscript o refers to that of an infinite cylinder 


= Reynolds number, D,Vp,/u 


thickness of hypothetical tube of gas 

temperature; subscripts have the following meanings: g the 
temperature of the ambient gas, p that of the particle, s the 
uniform temperature of the surface to which heat is trans- 
ferred, w the temperature of the furnace wall, © the uni- 
form temperature of the gas at infinite distance from the 
surface receiving heat 

lin = Urge 

absolute temperature, degrees Rankine; subscripts p and w 
refer to that of particle surface and wall surface, respectively 
components of fluid velocity in the x, y, z directions, re- 
spectively 

average volume of a particle 


= velocity of particles relative to furnace wall (in a stagnant 


gas); subscript ¢ refers to terminal velocity 


= average mass of a particle 


coordinates of a fluid element 


= V/ thermal diffusivity, = / k/cp 


volume coefficient of thermal expansion 
ratio of particle area to wall area 
temperature rise 


= absorptivity of cloud of particles, i.e., fraction of radiation 


incident on cloud that is absorbed 
emissivity of particle surface 


= roots of first order Bessel function 
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\ = characteristic linear dimension 

uu = viscosity of gas 

ax = ratio of circumference to diameter 

WY = conformal coordinates for potential flow around a sub- 
merged body of revolution. Note that VII is the velocity 
potential 

p = density; subscript p and g refer to the particle and the gas, 
respectively 

6 = time; Ad; = time of fall through furnace 

¢ = angular coordinate 

w = solid angle 


4. Mechanism of Heat Transfer —tThe transfer of heat from a solid 
surface through a gas to falling particles takes place both by con- 
vection and by radiation. The thermal resistance to convection is 
composed of the resistance from the wall to the gas and that from the 
gas to the particles, acting in series. Transfer of radiant energy takes 
place directly from the wall to the particles. If the gas has infra-red 
absorption bands, it also will receive heat directly from the wall and, 
consequently, it will partially blanket the particles from radiant 
energy. The net radiation transfer acts in parallel with the convec- 
tion. The thermal resistance within the particle is important only if 
the conductivity of the solid is low, or if the diameter is large. Con- 
fining our attention to the opposite conditions for the present, the 
overall coefficient of heat transfer is stated mathematically in terms of 
the individual coefficients by the equation 


Rep 
oe == h, e (1) 


hie 
ery 


yCW 


This equation apples only to point conditions. Each of the individual 
coefficients and the area ratio, y, varies with the position of the particle. 
In the case of particles falling through a cylindrical tube, y may be 
expressed conveniently in terms of a dimensionless group which also 
enters into the radiation coefficient 


aD? 
(c iL) (4a) 
NA, 4 


as = = (Da.* 2 
Ale arDdbL ‘ ( ) 


*The surface area of particles having convex surfaces is four times the projected area, cf. p. 32. 
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Little is known about the individual coefficients, as this problem 
in heat transfer has apparently received little attention. A discussion 
of the theory underlying the mechanism of heat transfer to falling 
particles will prove helpful in interpreting the experimental data and 
in extending the work to the practical design of large scale equipment. 


5. Heat Transfer to Small Particles by Convection—An exact 
analysis of this problem in convection is exceedingly difficult, if not 
impossible, because of the complicated hydrodynamics involved. Any 
solution of the problem must be made on the basis of simplifying as- 
sumptions which are valid at best only within certain limits. The re- 
sulting equations, therefore, must be recognized as approximations, 
and the nature of the assumptions should be clearly understood in 
order to avoid faulty application. 5 


A. Simple Concept 


One of the simplest concepts of the mechanism of heat transfer 
to falling particles is analogous to that recently proposed* for the 
absorption of gases by spray droplets. The assumptions, while ad- 
mittedly crude, lead to an expression for the heat transfer coefficient of 
the same form as that given by the more rigorous attack on the 
problem to be described later. Furthermore, the numerical values 
calculated in this manner agree fairly well with those obtained from 
the experimental data, so that, as a first approximation, the equation 
has been recommended for design purposes.f The corresponding 
equation for mass transfer was likewise found to agree within ten 
percent with the observed values of the coefficient for the absorption 
of ammonia, sulphur dioxide, hydrogen sulphide, and carbon dioxide 
by falling droplets of various solutions. 

In this concept, it is assumed that each particle makes contact with 
a tube of gas described by the trace of its periphery and of thickness s. 
This thickness represents the distance through which all of the heat 
above the temperature of the particle surface can be transmitted by 
conduction while the particular layer is in contact with the particle 
surface. As the velocity of the particle increases, s must decrease. In 
any case, s is small compared to the diameter of the particle, so that 
the volume of the tube is 7D,Ls. In its simplicity, the tube of gas 
represents a zone in which the temperature is reduced to that of the 
surface of the particle, since complete temperature equilibrium is 

*H. F. Johnstone and R. V. Kleinschmidt, Trans. Am. Inst. Chem. Engrs., vol. 34, p. 181, 1938. 


tUniv. of Ill. Eng. Exp. Sta. Bul. 324, p. 79, 1940. : 
tH. F. Johnstone and G. C. Williams, Ind. Eng. Chem., vol. 31, p. 993, 1939. 
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x 


established within the allotted time. The layex obviously does not cor- 
respond to the so-called “laminar film,” or “boundary layer” and is 
purely hypothetical. 

The quantity of heat transferred across the tube from the ambient 
eas to the particle surface as the particle moves through a distance 
of one diameter is 

2 
cae : (to ae ty) ae 0 
V 


From a heat balance for this segment of the tube, 


k p 
0 Os Gs = i) —= aD ,’spoC Ga cs bp) (3) 
s 
h aR 
whence ED; 
$ = (4) 
V pc 
The heat transfer coefficient h-., is 
k KV pgc 
in = ——-., (5 
7 Ss 4 D, ) 


With the exception of the proportionality constant, here equal to 
unity, the term on the right is indeed identical with that obtained by 
Boussinesq in his classical treatment of thermal conduction to bodies 
submerged in moving fluids. A brief discussion of the fundamental 
equation and of the customary assumptions required for its solution 
will now be given, along with a solution which avoids the most doubt- 
ful approximation in the Boussinesq treatment. 


b. Solution of the Fourier-Poisson Equation for a Spherical Particle 


The Fourier-Poisson equation is a differential equation for thermal 
conduction in moving fluids. In rectangular coordinates its form is 


at ot ot ot ) ot 0 at 0 ot 
Cp ie ) = (i ) + (i ) 4 (i \ (6) 
06 Ox Oy Oz Ox \ Ox Oy\ oy Oe OF 


This equation itself is derived from the laws of hydrodynamics and 
thermodynamics. In its simple form, as written,” it applies strictly 


site assumptions underlying the Fourier-Poisson equation and its solution for various cases 
are discussed by T. B. Drew, Trans. Am. Inst. Chem. Engrs., vol. 26, p. 26, 1931. 
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only to an ideal incompressible fluid in which the forces of viscosity are 
negligible, and in which the kinetie energy changes are small compared 
to the total heat change. These restrictions themselves do not seriously 
limit the usefulness of the equation, as for many fluids, either gaseous 
or liquid, moving through a heat exchange system at relatively low 
velocities, the forces of viscosity are negligible and the pressure is 
nearly constant. 

Application of Equation (6) to convection problems is made after 
determining, or assuming, the velocity distribution of the fluid, so that 
u, v, and w may be evaluated at all points in the system. The simplest 
case is that in which the velocity in all directions is zero. Such a con- 
dition would prevail approximately in heat transfer from a quiescent 
fluid to a very small particle. For the steady state in this case the 
solution®* of Equation (6) requires that 


f= lah Dat ee 3.) (7) 


Thus, the heat transfer coefficient is 


2k 
Wey = ap. (8) 


This agrees with the recent observation made from experimental data 
by Meyert that, for natural convection, the Nusselt number, h.,r/k, 
approaches a constant value between 0.5 and 1.0 for small values of 
the Grashof group, r?p?BgAct/u?. Meyer estimates that the heat trans- 
fer coefficient for a particle 0.0016 inch in diameter surrounded by 
air in a spray drier would be 240 B.t.u./(hr.) (sq. ft.) (deg. F.). With 
a temperature gradient of 200 deg. F., which customarily prevails in 
this apparatus, the rate of heat transfer would be 50 000 B.t.u./(hr.) 
(sq. ft. of particle surface). To obtain a rate of this magnitude by 
radiation alone would require a radiator temperature of about 2000 
deg. F. This emphasizes the high rate of heat transfer attainable for 
small particles. 
a. The Boussinesq Solution 


In order to develop the Boussinesq solution of Equation (6) for 
spherical particles submerged in a moving fluid, it is desirable to con- 
sider first the case of the submerged flat plate. From this the 


*The solution is best accomplished by writing the equation in spherical coordinates. 
+P. Meyer, Trans. Inst. Chem. Engrs. London, vol. 15, p. 127, 1937. 
tJ. Boussinesq, J. math. pures et appl., vol. 1, p. 285, 1905. 
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more complicated problem can be attacked“ by making additional 
assumptions. 

Let the flat plate of width x, and of infinite length be located in 
the x, y plane with the leading edge at the y axis. The fluid moves in 
the x direction parallel to the surface. It is supposed that v and w are 
everywhere and at all times equal to zero, that uw is constant and in- 
dependent of z, and that k, c, and p are constants. These conditions 
limit the case to streamline flow over the plate, although the results 
should be valid as an approximation even when there may be turbu- 
lence near the trailing edge of the plate, but when the boundary layer 
as a whole is essentially laminar. 

For the steady state, the foregoing conditions reduce the left side 
of Equation (6) to the single term, cpudt/dx. Boussinesq simplified 
the equation further by assuming that the second partials with respect 
to x and v are negligible. For a long plate it is obvious that the flow 
of heat in the v direction must be negligible and, therefore, 0?t/dy? 
will vanish. Neglect of the term 0?t/dx?, however, requires that the 
thermal conduction in the direction of flow also be negligible. This 
assumption will be shown to hold only within certain hmits of the 
fluid flow defined by the Reynolds number. 


The simplified equation now becomes 


ot a —O% 


ae . 9 
Ox u 02" © 


where a? is the constant k/ep. 
A solution of Equation (9) which satisfies the boundary condi- 


tions, t = t,, for all positive values of x and for z = o, and t = ¢,* 
for z = O for all values of x along the surface of the plate, is 
= tes 1 2 i ed 
2 =1- oP dk (10) 
iy = thes ~/ 1 0 
where 
P ae Neat 
2N az 


The rate at which heat is transferred from the fluid to the plate 


_ *Although the equation may be solved for any variation of ¢ with z, only the simple case 
of constant surface temperature need be considered here. 
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is found from the value of the temperature gradient at the surface. 
From Equation (10) 


ot 2 U 
Gye ee ay 


For the differential area dxdy located at the distance x from the 
leading edge 


ot U 
dq = —k (=) CRO = (te — 1) dxdy. (12) 
@ 2—0 


Torr 


The overall rate of heat transfer to the strip of width dy and length 
21 is 


“A. U 
C= 0 (tate) ay f _ : dx (13) 
0 CEi6b 


kcpu 
2k (ts — to) xidy ——., 


TX 


The average coefficient of heat transfer over the strip is 


2 kcepu kepu 
ie S mnee 2 Digs. 113, aa 
(ts — too)aidy Va N Ly 


Before proceeding further, it is well to examine the simplifying 
assumption made by Boussinesq. The assumption of negligible 
conduction in the direction of flow is valid only when the ratio of 
a*t/dy? to 0*t/d2? is large. 

If these derivatives are evaluated from Equation (10), the ratio 
of the first to the second is found to be 


oS) 7 
@ 3 ( eit ) Li ) 
G . ‘of — . 
2 


Thus the approximation becomes more exact the larger the values of 
Re- Pr and of x/z. It is not necessary to consider small values of 
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z/z, however, because most of the heat tran8fer takes place near the 
plate. If only that region is included for consideration in which the 
temperature gradient d¢/dz is greater than one-tenth of its value at 
the surface, and if the smaller of the two derivatives may be neglected 
when its absolute value is one-tenth that of the larger, the approxima- 
tion is valid for values of Re- Pr greater than eight. The principal 
assumption on which the Boussinesq solution is based is justified for 
diatomic gases when the Reynolds group exceeds 10.8, and for water 
(at ordinary temperatures) when the group exceeds 1.2. 

Extension of the Boussinesq solution for Equation (6) to a cylinder 
submerged in a moving fluid is readily made by transformation from 
rectangular coordinates to the orthogonal curvilinear system repre- 
sented by the stream function and the velocity potential of classical 
hydrodynamics. The only additional assumption required is that the 
motion of the fluid be irrotational, which in itself indicates that a 
velocity potential exists. By a suitable choice of the origin of coordi- 
nates, and after making the Boussinesq assumption of negligible con- 
duction along a stream line, the form of the equation becomes identi- 
cal with Equation (9). The boundary conditions are likewise the 
same and the solution is, therefore, similar to that already given for 
the flat plate. 

For the three-dimensional case of a submerged body of revolution, 
the form of Equation (6) in the new coordinate system now described 
by stream lines (¥ = constant) and equipotential surfaces (VI — con- 
stant) is 


ot a ot 2 oll ot Ot 
== = ae | (15) 


oll? (=) ) (=) dx OW ow? 
Ov Oy 


The x axis is the axis of revolution, and is parallel to the direction of 
flow, while y is the radial distance from the x axis. 

In order to simplify Equation (15), the first and second terms in 
the bracket are assumed to vanish and the value of y is limited to its 
value at the surface of the body. Neglecting 02¢/d1? represents the 
usual assumption of. negligible conduction along the stream lines, but 
the neglect of the second term and the limitation imposed on the value 
of y introduce new assumptions which must be examined. Solution of 
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the simplified equation as before gives the average heat transfer co- 
efficient for the general case 


KkVpc 
hee = u| = (16) 


where M is a constant dependent on the shape of the surface and A is 
a characteristic linear dimension. For the sphere, the average co- 


efficient is 
2 kV pc kVopc 
hae a ys 1134) ae (17) 
V4 D; Dy 


We note that the coefficient is slightly higher than that calculated 
from the less rigorous Equation (5). As a matter of fact, the experi- 
mental data obtained in this study actually indicate that the propor- 
tionality constant should be less than one. Furthermore, Lévéque’s* 
critical examination of the simplifying assumptions mentioned in the 
foregoing throws considerable doubt on the validity of the generalized 
Equation (16) for the transfer coefficient. It appears that the second 
term in the bracket of Equation (15) is actually negligible over most 
of the region surrounding the body in which the temperature gradient 
exists. On the other hand, the limitation on the value of y becomes 
an allowable approximation only when 


Il, — Ip) & 
Yar > 2) (18) 


Here Il, and Ip are, respectively, the values of the velocity potential 
at the trailing and leading edge of the body. For potential flow 
around a sphere it can be shown? that 


SDs 
= (19) 
2 


I — Ib = 


Also, for a sphere, Yu» is approximately D,/z. Substitution of these 


*A, Lévéque, Ann. mines, vol. 18, Ser. 2, p. 249, 1928. mAs F 
tL. Prandtl and O. G. Tietjens, ‘“‘Fundamentals of Hydro- and Aeromechanies,” p. 161, 
McGraw-Hill Book Co., New York, 1934. 


¥ 
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Fig. 1. GrometrricAL REPRESENTATION OF PROBLEM OF 
Hear TRANSFER TO A SPHERE 


values in the inequality shows that the Boussinesq assumptions are 


allowable when 
DEV IDY AW. c 
e =( =) ( ©) > 13 400. (20) 


a MK 1M 


In such ranges of flow there is little doubt that the entire analysis will 
break down due to the departure from ideal flow conditions. 


b. Solution in Spherical Coordinates in Terms of Bessel Functions 


The Boussinesq solution of the Fourier-Poisson equation for heat 
transfer to a sphere assumed that a velocity potential exists for the 
flow around the body. By making a slightly different assumption in 
regard to the flow pattern it is possible to solve the equation directly 
in terms of Bessel functions which can be evaluated. 

For the steady state, Equation (6) may be written in the spherical 
coordinates, 7, 6, and g (see Fig. 1) as follows: 


Oh - Sek OP eect 


or ip (ol) rsin@ do 


lo 6 ot 1 a) ; ot 1 07t 
G a leet ate == isin ) (21) 
On or r> sin @ 06 00 je Say) Neck? 


See 
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vr, Ve, and vp are components of velocity in the direction of increasing 
r, 0, and g, respectively. 

We shall assume that the velocity of the fluid around the sphere 
is everywhere parallel to the surface and equal to V. As a result of 
symmetry about the axis of revolution, v,;, vz, and 0°t/d¢g are zero, 
and ve = V. The second term in the bracket represents the ac- 
cumulation of heat in an element of fluid due to heat transfer in 
the direction of increasing 6. Since the streamlines have been as- 
sumed to be parallel to the surface, this represents the net conduc- 
tion along a streamline, and should be negligible according to the 
usual Boussinesq assumption. The equation then becomes 


a & 9 at. 
= (= i (22) 
00 Viarwor or 


A solution of Equation (22) may be assumed to have the form 


ee ee) (23) 


too — te 


where F (6) is a function of 6 alone and G(r) is a function of r 
alone. Differentiation of Equation (23) and substitution in Equa- 
tion (22) gives 


(24) 


anal a E dG | CG | 

r ‘ 
HY ao VG dr dr? 
Since each side of this equation is a function of only one independent 
variable, each may be set equal to a constant, say —b*. Then the 
forms of the functions F’(@) and G(r) are found by solving the ordinary 
differential equations: 

dF 


+O?Fr =0 (25) 
dé 
and 
2 IG b2V 
Beg Gat: (26) 
dr? dr a 


The first of these has the solution 


F(@) = Be”? (27) 
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XY 


and the second, ob 
bv/ rV a 


where J; is the Bessel function of the first kind and of the first order. 
Since the original equation is linear in G and its derivatives, the 
sum of any number of particular solutions is also a solution, or 


as = > Be,” 


et ee biv/ rv ci 


— a z — rV ) 


a 


The arbitrary constants B; and b; may be determined from the 
boundary conditions which require, first, that ¢ = ¢, for 6 > 0 and 
r =r, and, second, that.t =%, for @= 0 and fr 71. Prom the 


first condition, 
2b; rV 
eh I (30) 


a 
Thus, the successive values of the constant become 


nv N2 
; ; 5 EUG: 
IRQ 2 Ip VIRB 0 IPP 


ma 

acy 
I 

S 


The values of », which are the roots of Ji(¢) = 0, are given in 
standard tables.* 


From the second boundary condition 


(31) 


a 


2bivV PW 
Ca ee 1 d ) 
i=0 bia/ rv. rV 


The values of B; may now be determined by comparing the coeffi- 
cients with those of similar terms of a Fourier-Bessel seriest for the 
expansion of ~/ 7: 


Nes a (a8 n( zs a ; (32) 


Qa 


pie Jahnke and F. Emde, ‘Funktionentafeln,”’ Penbnery eee 1909. 


+F. Bowman “Introduction to Bessel Functions,” p. 8-11, Longmans, Green and GCo., 
London, 1938. 
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Fig. 2. Comparison or Equations or Hear TRANSFER TO A SPHERICAL PARTICLE 


The final solution of Equation (22) then becomes 


t—t, oo a eb: 2b rV 
pee ey yy a Sy | 
lhe =p i=0 Q2bir/ rV V/ rV 

bid { ——————_ 


a 


. (33) 


The total rate of heat transfer across the surface of the hemi- 


sphere is 
T T at 
q= -{ i k (—) r? sin 6 dgedé 
0 0 or TT, 


o jot enti 
aie See a 34) 
ee, Peay 


I 


The average heat transfer coefficient then is 


k 2 Lem? 


hep = (35) 
-o be x eee 


The series converges fairly rapidly. For Re-Pr=—10 the sixth term 
is only 0.1 per cent of the sum, while for Ae:-Pr— 100, seventeen 
terms are required for this percentage. 
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Equation (35) is plotted in Fig. 2 along with Equations (17) and 
(36). It is interesting to note that, according to Equation (35), the 
Nusselt group, hepD,/k, approaches the value of 2.0 as the Reynolds 
number approaches zero, as required by Equation (8) for the heat con- 
duction to very small spheres in a quiescent fluid. At high Reynolds 
numbers the curve approaches the value given by 


kV pe 
Rep = ond jo (36) 


Pp 


This is a limiting value which may be derived by integrating Equation 
(13) for heat transfer to a submerged flat plate over a spherical sur- 
face. Such an integration is justifiable when the boundary layer is 
very thin. In such a case the surface of the sphere may be conceived 
of as being made up of many narrow flat strips each of angular width 
df, and extending from the leading point to the trailing point. The 
diagrammatic representation of this problem is also shown in Fig. 1. 
We again assume that the velocity of the fluid elements is everywhere 
constant and parallel to the surface of the sphere. From Equation (14) 
the rate of heat transfer to an element of area, 7,? sin 6d6d¢, located 
at a distance, 7,4 from the leading point of the sphere is 


( = J cepVk 7, sin 6 déde (38) 
ag = Ga Uso = O 
: T 4/ r0 


Integration with respect to 6 may be performed after expansion of 
sin # in a series. This gives 


on Ee Wee) cpVk 
A/a NED 


q = (t, Qa? (39) 


and the average heat transfer coefficient becomes that given by 
Equation (36). 

Not only do the values of h-,» from Equation (35) agree’ more 
nearly with the experimental results than the Boussinesq Equation 
(17), but also there is reason to believe that the assumed pattern of 
flow around the sphere corresponds more nearly to actual conditions 
in the region near the surface in which the resistance to heat transfer 
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Fig. 3. ASSUMED CoMPONENTS oF FLUID VELOCITY IN 
NEIGHBORHOOD OF A SPHERE 


exists* than that required by the theory of potential flow. The com- 
ponents of velocity are compared with those calculated for potential 
flow and for viscous flow around a sphere in Fig. 3. The latter are 
derived from the Stokes’ equations; which neglect entirely the forces 
of inertia. At a point on the surface 90 degrees from the leading point, 
according to the potential flow theory, the fluid has a tangential 
velocity of 1.5V, while, according to the theory of viscous flow, the 


*The effective thickness of a uniform fluid layer surrounding a sphere in which all of the 
resistance to heat transfer may be considered to exist, expressed as a fraction of the radius of 
the sphere, is approximately 2k/hep Dp. According to Equation (35) this group is never greater 
than one and at a value of Re:Pr equal to 100 the value is 0.3. Therefore the distance in 
question cannot be greater than the radius of the sphere. 

tLeigh Page, “Introduction to Theoretical Physies,” p. 233, D. van Nostrand, Inc., New 
York, 1928. 
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velocity must be zero. Actually the velocity is “ero at the surface and 
increases rapidly within a narrow region to a value nearly equal to V. 
The radial component is zero according to either theory. Thus, at this 
location, the velocity has approximately the value assumed for the 
foregoing solution. At other points on the surface the assumption of 
the constant tangential velocity is not far in error. 


6. Radiation to Clouds of Particles—Heat transfer by radiation 
from a solid surface to small opaque particles is similar in principle to 
radiation to an absorbing gas. When a radiant beam enters a cloud of 
particles, those particles near the front of the cloud have a greater 
chance of intercepting the beam than have those in the interior or at the 
rear of the cloud. This is equivalent to saying that the average inten- 
sity of the radiation varies through the cloud. It is necessary to know 
the variation of intensity with distance along the beam in order to cal- 
culate the transfer of radiant energy to particle clouds of various 
shapes. The problem is analogous to that arising in gas radiation. In 
fact, the equations used for the calculation, with the introduction of 
new constants, will be similar to those used by Nusselt* for calculating 
radiation to different gas shapes. 

For particles which absorb the incident radiation completely the 
problem of variation of intensity with distance is that of finding at 
any point the average fraction of the beam blanked out by particles 
along its path. This is equivalent to finding the covering power of a 
number of small dises tossed at random on a plate. The analysis has 
been made by Haslam and Hottel} as follows: 

When n dises, each of area a, have been tossed on a plate of area 
A, the fraction covered will be f,. The probability that the (mn + 1)st 
dise will fall on uncovered area is 1 — f,. On the average, its effective 
covering power will be a (1—f,). For n+ 1 dises the fraction 
covered is 


A 


Z 


oe fn) (= jie 3 
n+1 — Jn n <) = tef Sate (39) 


In a similar manner, the fraction covered by (n + 2) dises is 


a a ay Ne a ONG 
= | 
Fnte = fat (1 ) = f.(1 - ; ) ee 
A A A A A 
*W. Nusselt, Forschungsarbeiten Geb. Ing., No. 264, p. 40, 1923: | Ter. De _ 
70, p. 763, 1926. g € ng ( i ; 3; Zeit. Ver. Deut. Ing., vol. 


tR. T. Haslam and H. C. Hottel, Trans. Am. Soc. Mech. Eng. (FSP 50-3), vol. 50, p. 9, 1928. 
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or 


fon 1=f(1-—) -[1- + (LY ]=u.-n(1-). 0) 


By mathematical induction, for (m + n) dises 


fain — 1 = (f,—) (1-2). (41) 


eee € S =). (42) 


If C is the number of solid particles per unit volume in the space 
above A and | is the thickness of the cloud above A, the number of 
particles, m, projected on A is CIA. Therefore 


ea (1 ~)o (43) 
Ue A ° 


When n = 0, 


When the area of each particle is very small compared to the area 
of the radiating surface, 


NC. 
ri -(0—4)"J 
a AlaqCla 
=l]-—- | Lim (2 — =) | 
a/A— 0 A 


= 1 — ee, (44) 


From Equation (44) the relationship between the average intensity 
at point / and the intensity at / = 0 is 


T, => T,e—¢'. (45) 


For particles which do not absorb completely, but which have a 
high emissivity so that the second order radiation is negligible, 
Equation (45) becomes 

Lopes deena, (46) 
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Fig. 4. GromerricaL ANALYSIS OF RADIATION FROM WALLS OF AN 
INFINITE CYLINDER TO CLOUD OF PARTICLES 


This equation is similar to Beer’s law for the absorption of light 
by liquids and gases except that here the proportionality factor may 
be calculated from the size of the particle, while in Beer’s law it must 
be evaluated experimentally. 

Equation (46) may be used for calculating the absorptivity of a 
cloud of particles of a given shape and concentration. Referring to 
Fig. 4, the rate at which heat radiating from a surface element dA, 
strikes an element of volume dV, located in a cloud of particles at a 
distance / from dA and subtending the solid angle dw at dA is given by 


dq, = I, cos a e~ "dA dw. (47) 


The average intensity of the ray at the front face of the volume ele- 
ment is obtained by dividing Equation (47) by the cross sectional 
area (l?dw) of the ray at the point l. 


COS @ Ee  &rCla 


han ee (48) 
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Now consider that the element of volume, dV, is subdivided into 
cylindrical elementary rays, each having a length dl and a cross sec- 
tional area dA’. The quantity of radiation entering each of these ele- 
mentary cylinders has the value Z;dA’. From Equation (46) the frac- 
tion ¢,Cadl of this amount is absorbed on passage through the cylinder. 
The total amount dq. absorbed by the volume element dV is eal pay. 
or 


€pCaln cos a- e~""dAdV 
dq. = p . (a2) 


The total radiation from a solid wall which is absorbed by a cloud 
of particles is obtained by integration of Equation (49) over the area 
of the radiation wall and the volume of the cloud, 


cos a: e~&l4d@ Ad V 
Op = eCal. | P ; (50) 


The total rate at which energy is emitted in all directions from the 
surface is I,7A. The absorptivity of the cloud of particles is the frac- 
tion of this total emission which is absorbed, that is, 


(! cos e-€ pCla 
& = “Sy = ae / Se 
nv area 


The integration of Equation (51) has been performed by Nusselt* 
for the case of the sphere and the case of the infinite cylinder. For the 


former, 


i 1 + 2€,Croa 
anes ap Po eter, (52) 
2 (€,Croa)? 2 (e,Croa)? 


For the infinite cylinder 
—2e,Cr a cos ¥ 
= l|- =" Als: — cs@ cos? 8 cosy dBdy. (53) 


The integral in Equation (53) was evaluated graphically. 


*W. Nusselt, Zeit. Ver. Deut. Ing. vol. 70, p. 763, 1926. 
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TABLE 1 Ze 


ABSORPTIVITY OF CLOUDS OF PARTICLES IN FURNACES oF VARIOUS SHAPES 


Absorptivity €c Absorptivity €¢ 
2epCroa Between 2epCroa ‘ Between 
; Sph Infinite Infinite nhere Infinite Infinite 
Re Cylinder Parallel P Cylinder Parallel 
Plates* Plates* 
0 0) 0 0) 3.50 0.8547 0.9293 0.9901 
0.10 0.0642 0.0932 0.1674 4.00 0.8865 0.9460 0.9945 
0.20 0.1243 0.1767 0.2961 4.50 0.9073 0.9580 0.9969 
0.25 0.1520 0.2149 0.3506 5.00 0.9232 0.9665 0.9983 
0.30 0.1792 0.2512 0.3999 5.50 0.9356 0.9728 0.9989 
0.40 0.2308 0.3170 0.4854 6.00 0.9454 0.9776 0.9994 
0.50 0.2748 0.3764 0.5568 6.50 0.9532 0.9814 0.9996 
0.75 0.3835 0.5000 0.6905 7.00 0.9595 0.9843 0.9998 
1.00 0.4715 0.5957 0.7806 7.50 0.9646 0.9868 0.9999 
1.25 0.5451 0.6709 0.8428 8.00 0.9688 0.9885 1.0000 
1.50 0.6066 0.7298 0.8866 8.50 0.9724 0.9901 1.0000 
ee 0.6596 0.7769 0.9175 9.00 0.9753 0.9913 1.0000 
2.00 0.7030 0.8142 0.9397 9.50 0.9778 0.9927 1.0000 
2.50 0.7719 0.8689 0.9674 10.00 0.9800 0.9940 1.0000 
3.00 0.8220 0.9048 0.9822 co 1.0000 1.0000 1.0000 


*Here ro represents the half distance between the plates. 


For the case of infinite parallel plates, Jakob* obtained the equation 


a/2 —2€,Crya 
Ee = 1-— 2 | GW ge Sia eo COS a: da. (54) 
0 


Here 7, represents one-half of the distance between the plates. The 
integral in Equation (55) may be evaluated in terms of the exponential 


integral 
YD e-z 
Ei(—p) = | dx 
0 4 0) 


for which values are listed in standard tables.; The equation then 
becomes 


éo = 1 — (1 — QepCroa) e222 + (2e,Croa)?Ht (—QepCroa). (55) 


All of these shapes are of practical importance in furnace design. 
The values of «. for the three cases are listed in Table 1 and are plotted 
in Fig. 5. Inspection of the curves shows that for the same value of 


*M. Jakob, “Der Chemie-Ingenieur.’”’ vol. 1, No. 1, pp. 299-303, Akad. Ve >1pzi 
tJahnke-Emde, “Functionentafeln,” loc. eo rh Bee Eee 
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Fic. 5. ABsorPTiviry or CLoups or ParticLes IN FURNACES OF VARIOUS SHAPES 


2e6,Croa, or e,CDa, a cloud of particles enclosed between infinite 
parallel planes has the greatest absorptivity. This may be explained 
qualitatively by the fact that the average length of a radiant beam is 
longest for the case of infinite parallel planes because of infinite di- 
mensions in two directions. The absorptivity of the cloud enclosed 
in an infinite cylinder is higher than that for a sphere because the 
cylinder has infinite dimensions in one direction. The absorptivity for 
other shapes of clouds encountered in practice may be estimated from 
Fig. 5 by interpolating between the curves shown. For example, the 
value for a cube should fall between the curve for an infinite cylinder 
and that for a sphere, since the average length of a radiant beam is 
less than that in the cylinder but greater than that in a sphere. The 
absorptivity for a long shaft of rectangular cross section should fall 
between the two upper curves.” 


*In a recent paper, Hottel and Egbert show that the effect of the gas shape on radiation 
of furnace gases may be approximated by using an effective average beam length. On this basis, 
the absorptivity of a hemisphere of a cloud of particles receiving radiation from a spot on the 
center of the base of the hemisphere of radius 7 is 1— e-€rCre, When the gas shape is other 
than a hemisphere, the same formulation suffices provided the radius of the hemisphere is 
replaced by the effective beam length, 7, which, for small values of the quantity epCra, is four 
times the mean hydraulic radius, 7.e., four times the gas volume divided by the area of the 
bounding walls. This is exact as €e approaches zero. Comparison with the theoretical curves of 
Figure 5 shows that the suggested approximation is good up to values of €¢ of 0.2 but above 
this value considerable deviation occurs. For high values of the group analogous to €pCra in 
gas radiation, Hottel and Egbert recommend using 85% of the effective average beam length 
obtained as above. This appears to be satisfactory for €c greater than 0.6. While this method 
of approximation has the advantage of conforming to established engineering practice in the 
calculation of heat transfer from radiating gas masses, it lacks the theoretical background of 
Equations (52), (53) and (54). (H. C. Hottel and R. B. Egbert, Preprint of paper presented 
before Am. Soc. Mech. Eng., New York, Dec. 1940.) 
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. ayers . . “ 
The area-emissivity factor 4, is given by 


1 
Ny (56) 
1 1 
Se enh 


Ew Ec 


where ¢, is the absorptivity of the wall. Since «- is generally small 
compared to €,, Faz may be taken as equal to «-. The coefficient of 
radiant heat transfer to clouds of particles expressed on the basis of 
the area of the particles then is given by 


0.172, |(T,,/100)4 — Ee LOO t 
1, = BT [(Po/100)¢ — (7,/100)4]_ - 
oY Ci = i) 


7. Heat Transfer from Wall to Gas—When a cloud of particles 
falls through a heated tube, three factors produce convection within 
the tube, viz., the normal convective forces due to density differences 
along the wall, the increased convection due to the introduction of cold 
particles at the top of the tube, and the disturbance due to the passage 
of the particles through the gas. Consequently, the rate of heat 
transfer from the wall to the gas is somewhat greater than that due to 
natural convection. Because of the complexity of the problem no at- 
tempt has been made to reduce it to a mathematical basis. 

In the summary of the work previously published,* the Rice equa- 
tiony for the heat transfer coefficient for natural convection from a 
long vertical pipe was recommended as a rough approximation. In 
this the coefficient he» is given as a function of the Grashof group. 


k [ D3p,298 (tw — ty) 7/3 
hee ORES =| po 9B ( 2] ; 


: (58) 
2 

It was realized that the values of he» from Equation (58) would be 
quite conservative. From the treatment of the experimental data 
which follows it appears that the actual values are from two and one- 
half to four times those calculated from the Rice equation. The co- 
efficient 1s apparently independent of the wall temperature and of the 
particle size but increases with the diameter of the furnace. 


*Univ. of Ill. Eng. Exp. Sta. Bul. 324, p. 74, 1940. 
+W. H. McAdams, “Heat Transmission,’ p. 254, McGraw-Hill Book Co., New York, 1933. 
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TABLE 2 
PrincipaAL Dimensions or Furnaces Usep In ExpERIMENTAL Work 


Small Furnace Large Furnace 
Insidesdiameterlofitubesims. occaacuiee naccs de eaecs 1.61 3.76 
@utsideidiameteriof tubes imiaac. 4000. acres sine css 1.90 4.00 
Inside diameter of refractory lining, in............... B25 5.0 
Motaldencthvol tulbeninawyen cet iae nies ate eect 54.5 120.0 
Heated: lenothtot tuberin: . iiss aiciecc ses slersetat ecm 48.0 110.0 
Distance from top of calorimeter to bottom of heated 
DECULOM eA Ds ticitela ren) cori eon atcus <5 eierc winch ete reed crate Mew ee 2.0 4.0 
Distance from feeder to top of heated section, in...... 4.0 17.0 
Location of thermocouples above bottom of tube, in... 8, 20, 32, and 44 27.6, 49.2, 70.8, 
+ 92.4, and 114.0 


III. ExperRImMeNTAL WorRK 


8. Introduction—The object of the experimental work was to 
provide accurate data covering a sufficient number of conditions so 
that the individual heat transfer coefficients could be determined. The 
agreement of these with the theoretical equations or with empirical 
correlations should justify extrapolation to large scale design. To 
this end, measurements were made on the rate of heat transfer to 
particles of known average sizes of several materials falling through 
heated tubes at known wall temperatures. The rate at which the 
particles were introduced was varied over a considerable range as this 
proved to be the important variable by which the data could be 
analyzed to calculate the individual coefficients. The effect of changing 
from air to carbon dioxide as the ambient gas was also studied. 

The general method of the study consisted in adjusting the walls 
of a tubular furnace to constant temperature and then introducing 
the fractionally screened particles at a constant rate. The quantity of 
heat transferred and the temperature of the particles were determined 
by collecting the material in a calorimeter immediately below the 
furnace and observing the temperature rise of the water therein. The 
particles studied were sand, carborundum, and aloxite. These materials 
were selected because they are inert to air at high temperatures and 
do not react with water. Furthermore, their heat conductivities and 
absorptivities differ considerably. 


9. Heat Transfer Measurements—Two sizes of furnaces were 
used. The essential dimensions of these are given in Table 2. Details 
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Fig. 6. Apparatus Usep ror Heart Transrer MBbASUREMENTS 


of the small furnace and the top of the large furnace are shown in 
Fig. 6. In both cases stainless steel tubes were used and thermocouples 
were peened into the walls from the outside. The smaller furnace was 
heated electrically by resistance coils wound in four sections on a 
corrugated alundum sleeve. Uniformity of temperature was main- 
tained by adjusting the external resistance of the heating elements. 
A typical temperature distribution along the inside of the tube is 
shown in Table 3. 

The large furnace was heated by natural gas. The stainless steel 
tube was surrounded by cylindrical fire-brick tile, 1 in. thick, for pro- 
tection against the direct action of the flames. The outer wall of the 
furnace was constructed of 4%-in. fire-brick with nine burner ports 
arranged for tangential firing. Control of this furnace was manual. 
The temperature of the tube wall indicated by the six thermocouples 
was recorded by a multipoint recording pyrometer. Variation along 
the tube did not exceed 25 deg. F. during a run. 

The feed to the small furnace was introduced into the center 
through a glass tube. The rate was controlled by an electric vibrating 
feeder. The feed to the large furnace was introduced from a hopper 
by means of a motor-driven screw. This discharged into a water- 
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TABLE 3 
TypicaAL TEMPERATURE DISTRIBUTION IN SMALL FURNACE 


Distance Below Top of Temperature Distance Below Top of Temperature 
Heating Zone, in. deg. F. Heating Zone, in. deg. F. 

2 707 26 1047 

5 992 29 1052 

8 1105 32 1072 

11 1130 35 1092 

14 1085 38 1112 

17 1075 41 1107 
20 1075 44 1030 
23 1067 47 812 


Average temperature, omitting top and bottom points = 1075 deg. F. 
Average temperature from four wall thermocouples = 1080 deg. F. 


cooled section of vertical pipe fitted in the ‘center of the stainless steel 
tube through asbestos packing. This arrangement was necessary in 
order to avoid heating the particles in the screw feed by conduction. 

The calorimeter used was a four-liter Dewar flask fitted with a 
wooden top through which a stirrer was attached. With both furnaces, 
the tube extended about one inch through a hole in the top of the 
calorimeter. A forty-junction copper-constantan thermocouple, with 
one set of junctions in ice, was used to measure the temperature rise. 
The thermocouple was carefully calibrated against a Beckmann ther- 
mometer. A Leeds and Northrup students’ type potentiometer with 
extended lower scale was used to measure the thermocouple e.m.f. 
The water equivalent of the calorimeter, thermocouple, and stirrer was 
determined electrically several times during the work. 

The general procedure adopted for the runs was as follows: A 
known small quantity of the powder was introduced into the hopper 
of the feeder. When the temperature of the furnace was uniform and 
constant and the rate of temperature rise of the calorimeter was con- 
stant, the feed was started. The rate was maintained constant until 
the hopper had discharged to a fixed mark. In order to avoid a draft 
action through the tube, the hopper was not allowed to empty com- 
pletely. The elapsed time was measured with a stop watch. The net 
increase in temperature of the calorimeter was noted after the rise 
had again reached a constant rate. The wall temperature was taken 
as the average of that observed just preceding and just after the run. 
In general, these did not differ by more than 5 deg. F. 

Measurements were made in this way in the small furnace at wall 
temperatures from 575 to 1050 deg. F., and at feed rates from 3 to 
50 Ibs. per min. per sq. ft. Three particle sizes of sand and two of 
carborundum and aloxite were used. Two series of runs were made 
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with sand at wall temperatures of 725 and ‘1050 deg. when the air 
in the furnace had been replaced by carBon dioxide. In the large 
furnace two series of runs were made at a wall temperature of 1000 
deg. F. using two screen fractions of sand. 


10. Secondary Measurements —In order to calculate the overall 
heat transfer coefficients from the measurements it was necessary to 
have accurate data on the average diameter, the projected area, and 
the rates of fall of the particles. This necessitated a secondary study 
in which these quantities were measured accurately by improved 
methods. 

A. Average Particle Size 


The arithmetic average “diameter” of the particles in each screen 
fraction was determined by observation under a petrographic micro- 
scope equipped with a movable stage. The “true average diameter” 
was obtained by securing random orientation on the microscope slide 
by using a method suggested by Tooley.* A thin layer of a gum arabic 
solution in glycerine and water was spread uniformly over the slide. 
Pin scratches in the layer were observed under the microscope until 
the film reached a stiff jelly-like consistency. The material to be 
measured was mixed thoroughly and a small quantity was taken from 
the whole with a spatula. This was spread on a paper which was then 
held a short distance above the slide and tapped lightly to transfer the 
particles with a uniform distribution. The orientation of the particles 
remained fixed as they fell on the gum. 

Several hundred particles on each slide were traversed under the 
microscope from several directions. The average diameter was found 
by dividing the reading on the micrometer screw by the number of 
particles measured in each traverse. The average for all the traverses 
was then taken as the “true average diameter,” D,. 


B. Average Projected Area 


Tooley* has shown that for particles having plane or convex sur- 
faces the surface area available for heat transfer is four times the 
average projected area obtained for all possible orientations of the 
particles. The slides prepared as described were placed on a projecting 
microscope and the images on a ground glass plate were traced on 
paper. The projected areas of individual particles were then measured 
with a planimeter. The areas of 50 to 100 particles from each sample 


*F. V. Tooley, Ph.D. Thesis in Engineering, University of Illinois, 1939. 
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TABLE 4 


Summary or MrAsurEMENTS or DiamMETER, Mass, AND Density or PARTICLES 
Usrep In Hear Transrer Work 


Average Projected Area a 
Average Average ; sq.in. X 1072 
Screen Diameter Particle Density oe 
Size Dp Mass wy gm. per ce. K 
mm. gm. X 104 Experi- Calculated 
mental wD y?/4 
A. Sand 
30-40 0.545 1.865 2.65 0.435 2.50 25 
40-50 0.436 0.941 setae 0.429 1.40 1.39 
60-80 0.326 0.35" Pee 0.431 0.904 0.904 
Av. 0.432 
B. Carborundum * 
No. 40 0.502 1.214 3.21 0.290 1.92 1.98 
No. 60 0.344 0.386 estas 0.351 1.02 0.935 
Av. 0.320 
C. Aloxite 
No. 46 0.485 1.362 3.92 0.304 1.88 1.85 
No. 60 0.344 0.581 sane 0.344 0.949 0.935 
Av. 0.324 


were averaged. This method of determining the surface area is 
obviously not applicable to particles having concave surfaces. With 
none of the materials studied was this condition serious. 


C. Average Mass and Volume of Particles 
The average mass of a particle was found by counting and weigh- 
ing from one thousand to three thousand particles out of each fraction. 
The average volume of the particles was calculated from the 
average mass and the density. The latter was found by displacement 


of water in a volumetric flask. 
From the average volume and diameter of the particle the “volume- 
shape factor” K was calculated from 


Up 


D,3 


K = (59) 


Values of the factor were approximately constant for different sizes 
of the three materials studied but were much larger for the sand than 
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for the carborundum or aloxite particles. In all cases the coefficient 


was lower than z/6 which would be the value for spheres. 
A summary of the data obtained in these measurements is given in 
Table 4. 


TABLE 5 


Summary or ResuLtts oN DraG ComFrricIpNTs FOR SAND, 
CaRBORUNDUM, AND ALOXITE PARTICLES 


Average Diameter | Observed Terminal Drag Reynolds 
Screen Size of Particles Dp Velocity Vz Coefficient Number* 
mm. ft. per sec. Cp DypVtpg/u 
A. Sand 
30-40 0.584 he ect NN eens 
0.534 12.4 0.933 128 
0.568 sare 0.941 139 
0.583 13.5 0.900 149 
0.594 14.6 0.875 162 
0.584 15.9 0.671 173 
40-60 0.401 St86° ilk tee © OL aes 
0.427 9.35 1.36 75.8 
0.435 9.88 1.25 81.5 
0.465 10.2 1.24 90.8 
0.490 10.8 es 100 
0.518 12.1 0.966 120 
60-80 0.313 i corel OO Me Senin Meme I omax 
0.308 6.05 2.54 34.9 
0.319 6.57 2.12 39.3 
0.335 7.08 1.92 44.5 
0.387 7.57 1.81 52.4 
B. Carborundum 
No. 40 0.464 SSN ere SOD) a ee eee 
0.490 9.00 1.52 79.8 
0.490 9.70 1.33 85.5 
0.483 Oe ile lye 89.1 
0.499 10.9 1.07 97.8 
0.523 pW es 0.988 109 
0.535 12.3 0.893 119 
No. 60 0.272 FisOD= 9 ie iets Ml Se aetna 
0.302 6.12 2.08 33.1 
0.326 We LS 1.75 40.5 
0.365 8.06 1.49 51.9 
0.381 8.79 Wid 59.9 
C. Aloxite 
! 
No. 46 0.480 TOR = |||) cern mI mune ore 
0.494 115 1.18 101 
0.507 125 1.02 113 
0.524 LS av 0.880 128 
7 0.591 15.7 0.791 162 
No. 60 0.331 EDO; Se © il) 9 - Riecepeeeee 9) meeellimae ners 
0.348 8.50 1.56 51.6 
0.371 9.58 1.31 62.5 
0.410 LOWE 1.16 Wve 
0.430 12.3 0.941 92.2 


*Based on physical properties of air at 85 deg. F. 
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D. Terminal Velocity and Drag Coefficient of Particles 


The terminal velocity of the particles was measured by air elutria- 
tion tests using a method similar to that of Martin.* The pro- 
cedure was as follows: About 25 grams of the material from each 
screen fraction was placed in a brass cone at the base of the elutriation 
tube. This was a glass tube 1 in. in diameter and 51 in. high. Air was 
blown into the cone through an inverted nozzle. The velocity was 
measured by means of a calibrated Venturi meter. Beginning at a low 
velocity the air velocity was gradually increased until a few particles 
were thrown a short distance above the top of the elutriator tube. At 
this velocity, the flow of air was maintained constant until only a very 
few particles were being ejected, and these were rising only a short 
distance above the tube. The air was then shut off and the material 
collected outside of the tube was weighed. An effort was made to 
adjust the air velocity so that each fraction had approximately the 
same weight. Several fractions were collected from each sample placed 
in the elutriator, and the average diameter of the particles in each 
fraction was determined as described in the foregoing. In calculating 
the drag coefficients, the arithmetic average of the air velocities of 
two successive sub-fractions was taken as the terminal velocity of the 
particle having the average diameter of the second sub-fraction. The 
coefficients were calculated from the equation 


SK Pp — Pog Dog 


79 


TT Pg et 


Cp = (60) 


It may be noted that in this equation the projected area of the particle 
on which the drag resistance is based is taken as 7D,°/4. There may 
be some question as to whether this is correct, since it has already 
been shown that the volume-shape factor was less than 7/6 for all of 
the materials studied. For convex surfaces, however, the error in- 
volved in computing the projected area as that of a sphere having the 
diameter D, is less than that in computing the volume on this basis. 
Furthermore, the projected area calculated in this way agrees closely 
with that measured directly, as is shown in Table 4. 

Values of the drag coefficient calculated for the various materials 
are shown in Table 5. The table also shows the variation in the size 
of the particles in each screen fraction used in the heat transfer 


*Geoffrey Martin, Trans. Inst. Chem. Engrs., London, vol, 4, pp. 164-78, 1926. 


Drag Coerficiers7, Cp 
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Fig. 7. Drag CorrricieENts For Partictes Fattinc THroucH AIR 


measurements. The drag coefficients are also plotted in Fig. 7. The 
data agree excellently with Martin’s results on carefully sized fractions 
of sand particles. When the volume-shape factors were not included, 
the points representing the data on the three materials were scattered 
considerably. Since inclusion of the factors brings the data together, 
the surface roughness evidently is not important in determining the 
frictional coefficient. The carborundum particles were considerably 
rougher than the sand, as shown by the smaller shape factor. Under 
the microscope also these particles appeared to have much more ir- 
regular shapes than the sand and aloxite particles. Martin’s volume- 
shape factor for sand was 0.284, this indicating very rough grains. 
It is probable, however, that he did not take precautions to prevent 
preferred orientation of the particles on the microscope slide and that 
the lower value may be explained on this basis. 

From the drag coefficients, the time of fall was calculated from the 
familiar equation of accelerated motion,* 


dV x One prey? 


Bese = 61 
ap: 9% rane ee ope ot 


*A slight error is recognized in this calculation since the effective mass of a body is greater 
for accelerated motion than for uniform motion due to the inertia of the ambient fluid. The 
mass is increased by approximately one-half of that of the displaced fluid. 
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TABLE 6 
CaucuLatep Times or Fatu or Parricyes Turoucu FuRNACES 


Time in Heated Section of Tube, sec. 
Screen Size Ambient Gas 
Small Furnace Large Furnace 
A. Sand 
30-40 Air 0.474 0.823 
40-60 Air 0.518 0.976 
40-60 COz (aes 7 i=" Weegee 
60-80 Air 0.613 1.258 


B. Carborundum 


No. 40 Air OV > i Saks 
No. 60 | Air ORGT ON ee | ne 
C. Aloxite 
No. 46 Air D835 | caw ae 
No. 60 Air O2D52— °° ~ 5 . sees 


Equation (61) was integrated graphically, or by a series approxima- 
tion in the range between Reynold’s number of 10 and 200 where 


Cp = 12.8Re-0-53, (62) 


The distance of fall as a function of velocity was hkewise determined 
by a second series integration. The time required for a particle to pass 
through the heated zone was then found from graphs of the velocity- 
time and distance-velocity functions. 

It may be observed that the product Cpg is nearly independent 
of temperature and, therefore, for any one gas, the time of fall caleu- 
lated for one temperature may be safely used at other temperatures 
not too far removed. The viscosity of air varies approximately as the 
0.7 power of the absolute temperature. Since the density varies as the 
~1.0 power, the Reynolds number varies as the —-1.7, and the drag co- 
efficient as the 0.9 power. Thus Cpp, must depend only on the —0.1 
power of the absolute temperature. In the experimental work, the 
maximum variation of the “film” temperature of the particles was 
from approximately 300 deg. F. to 800 deg. F. This would cause a 
variation in Cpp, of only 5.5 per cent. 

Table 6 shows the calculated average times of fall for the various 
materials studied in the two furnaces. 
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IV. CorRELATION OF EXPERIMENTAL DATA 


11. Experimental Results—The complete experimental data in 
the heat transfer measurements are given in Tables 7, 8, and 9. The 
runs are classified according to the furnace wall temperature, particle 
size, quantity of material fed, size of furnace, and nature of the 
ambient gas. It will be observed that the basic quantities measured 
by the calorimetric observations give the final particle temperature as 
well as the total heat transferred. The logarithmic mean temperature 
difference between the wall and the particle was calculated and, from 
this, the rate of heat transfer per degree difference was computed. 
On the basis of the feed rate, surface area of each particle, and the 
time of exposure, the overall heat transfer coefficient, designated as hy», 
was then calculated. The radiation coefficient, h,, was estimated on 
the basis of Equation (57). The overall convection coefficient, based 
on the particle area, is equal to (h», —h,). 


A. Sample Calculation 
A sample calculation of the data, based on Run 95, is shown below. 


Weighitvofiteed .a5-) sas Aste @ tere RN Ron ne ee ee 100 gm. 
Duration of feed = ..c.rut.s eo cene acme rote ona ee rene eee 80.3 sec. 
wl. ; T USL Nye 

Cross-sectional area of furnace.................. ( --) ( i ) = (0.0141 sq. ft. 

F, feed rate OS = 11.7 (Ibs.)/(min.) (sq. ft 
SLCCOLUALC pee apn imesh Baee0s C0 ol .7 (Ibs.) /(min.) (sq. ft.) 
/N, SUNG aba OP NHICLS HOTTA. 5 op aaokeeccanaonboungucescedoce: 361 deg. F. 
Meanrspeciticuheatrotoit aituzeiee si ra ae ee 0.205 
Heat transferred during rum, ..j..0..5...-++..40- ay X< 0:205 K 3.5) = 16:6 Bittu: 
Rateroneateunansic a ara yaa er ie = 745 B.t.u./hr. 

80.3 
isp) ENR UREN TENORS 5 ots ook dosahoonAoonudengsedcapasane 858 deg. F. 
(tw—tm), temperature difference at top of furnace...........858 — 75 = 783 deg. F. 
(tw —ty2), temperature difference at bottom of furnace. .858 — (75+361) = 422 deg. F. 
Ati, logarithmic mean temperature difference........... ee epee = 585 deg. F 
783 ee 
In : 
422 
; 6F Da AT 7 2.6 a 

yates talib, LOU DaXd 4 0.475 X 11.7 X 1.61 X 2.65 & 10-8 = polee 
60w,L 60-X 412 104 X48 oe2e 
Aw, Wal ATC a .is.ho da ow ade ase eee eee wDL = 1.686 sq. ft. 
NA,, total particle area in tube............ yAw = 1.686 X 0.0198 = 0.0334 sq. ft. 


16.6 X 3600 
80.3 X 585 X 0.0334 
= 38.2 B.t.u./(hr.) (sq. ft.) (deg. F.) 


Am, Overall heat transfer coefficient 
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EL ODSOLDLLVAb ya OLN UArh ze seen ern ane eee 0.5; (e,CDa) = 0.5 X 0.0198 = 0.0099 
<A SOLDUiVity OUcloud (rable: | yet er Bek iy esa acct fe. ee. 0.0099 
h,, radiation coefficient, based on arithmetic average particle temperature 


( 858 ony (2 + 460 \4 
0.172 * 0.0099 100 Ton ) | 


0.0198 585 


(hm — h,), net transfer coefficient due to conduction-convection exchange alone. . 
38.2 — 3.6 = 34.6 B.t.u./(hr.)(sq. ft.) (deg. F.) 


= 3.6 


B. Validity of Assumptions Made in Calculation of Data 


In the calculation of the heat transfer coefficient as described in the 
foregoing the following basic assumptions were made: (1) that a 
steady state was rapidly established between the wall and the gas, 
(2) that the average temperature of the particle as determined from 
the calorimeter measurements was equal to the surface temperature 
and (3) that the overall rate of heat transfer is proportional to the 
logarithmic mean temperature difference between the wall and the 
particle. These assumptions will now be investigated. 

The first assumption appears to be logical on account of the small 
heat capacity of the gas in the furnace and the high heat conductivity 
of the walls. As proof that a steady state between the wall and the 
gas was reached, several series of runs under various conditions were 
made at approximately constant rates of feed for varying lengths of 
time. These are shown in Table 7E to 7], inclusive. When the ob- 
served heat transfer coefficients are corrected for slight variations in 
the feed rate, the results indicate that the steady state was reached in 
less than ten seconds, as no significant variation in the overall co- 
efficient was found over periods from ten seconds to four minutes even 
at high rates of feed. 

The second assumption requires that the unsteady radial flow of 
heat within the particle be rapid compared with the rate of heat 
transmission up to the particle surface. The ratio of the thermal re- 
sistance outside the particle to that inside is equal to 2kh/hepD». Mc- 
Adams? shows that when this ratio is greater than about 6, the temper- 
ature gradient within the particle may be neglected. It will be shown 
later that the experimental values of the surface coefficient vary from 
about 50 B.t.u./(hr.) (sq. ft.) (deg. F.) for the largest particles to about 
30 for the smallest. Using the value 4.2 B.t.u./(hr.) (sq. ft.) (deg. F./ nue) 
for the thermal conductivity of quartz, which has the lowest value of 
any of the three materials studied, the values of the thermal resistance 


*K. Hild, Mitt. Kaiser Wilhelm Inst. Eisenforschg. (Dusseldorf), vol. 14, Part 200, pp. 


59-70, 1932. ; F : 
tw. H. McAdams, ‘‘Heat Transmission,’”’ p. 87, McGraw-Hill Book Co., New York, 1933 
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TABLE 7 % 
Heat TRANSFER TO CLOUDS OF SAND PARTICLES 


hm | hy | hin — hr 
F 
t) : tw A Alim NAp 
Run lb./(min.) = f 
sec. a ft.) deg. F. | deg. F. | deg. F. sq. ft. Be ays) Gaaitewanticle 
surface) (deg. F.) 
A. Small Furnace; Nominal Wall Temperature 1050 deg. F. 

a. Particle size—30—40 mesh; 100 gm. fed 
136 39.6 23.57 1074 414 GUS 0.0694 33.0 Selh PE 8) 
137 37.0 25.4 1072 410 775 0.0745 32.6 De 26.9 
138 Bone 26.6 1073 413 774 0.0780 33.0 bose Oss) 
139 38.0 24.7 1073 405 780 0.0724 32.0 Baa 26.3 
140 323 2.92 1066 530, 693 0.0856 48.5 6.2 42.3 
141 76.7 12.4 1058 444 739 0.0364 olen ‘sya h 31.4 
142 86.7 10.8 1042 428 733 0.0316 3653 SHO) 30.8 
143 210.8 4.45 1064 524 697 0.0130 48.0 6.1 41.9 
144 254.8 3.68 1057 SP» 690 0.0108 48.0 6.1 41.9 
145 125.8 7.46 1054 468 722 0.0218 40.7 5.8 34.9 

b. Particle size—40-60 mesh; 100 gm. fed 
126 25.4 37.9 1092 436 780 0.1330 28.2 5.4 2228 
127 40.4 23.2 1075 454 750 0.0835 30.8 6.0 24.8 
128 45.2 PADI 1069 463 738 0.0745 Soak 5.9 2652 
129 127.8 Tf so) 1072 570 672 0.0264 44.6 625. 38.1 
130 197.0 4.76 1067 600 6388 OLOLT 49.8 6.7 43.1 
131 108.6 8.65 1062 560 669 0.0312 43.8 6.3 Sh aD 
132 Done 17.6 1054 475 715 0.0634 34.1 Vf 28.4 
133 147.8 6.35 1045 580 636 0.0228 48.0 (Jac 41.7 
134 73.8 Pa 1044 509 684 0.0457 38.6 5.9 Sond 

c. Particle size—60-80 mesh; 100 gm. fed 
146 184.0 5.10 1051 686 583 0.038 29.8 6. 23.0 
147 76.4 1253 1045 593 618 0.092 24.8 6.3 18.5 
148 250.2 MD) 1039 678 558 0.028 31.8 6.9 24.9 
149 54.6 2 1027 558 633 0.128 22.3 5.9 16.4 
150 25.0 37.5 1067 461 738 0.280 15.4 Seats 9.7 
151 51.6 LEZ 1053 547 667 0.136 20.7 6.0 14.7 
152 50.4 18.6 1043 578 636 0.139 23.2 Gyal lyf pal 
1538 221.4 4,24 1033 639 580 0.082 28.4 6.6 21.8 
154 1225 7.66 1027 616 591 0.057 26.4 6.3 20.1 

B. Small Furnace; Nominal Wall Temperature 850 deg. F. 

a. Particle size—30—40 mesh; 100 gm. fed 
94 34.4 Dine 865 315 620 0.0796 30.2 3.9 26.3 
95 80.3 aN ay 858 361 585 0.0343 38.2 3.6 34.6 
96 148.0 6.34 852 385 563 0.0186 41.7 4.0 Otae 
97 128.5 7.30 849 371 569 0.0214 40.8 4.0 36.8 
98 226.9 4,13 846 420 533 0.0121 48.8 4.2 44.6 
99 140.0 6.70 833 376 549 0.0196 41.8 3.9 37.9 
100 86.0 10.9 827 371 545 0.0320 41.3 4.1 37.2 
101 100.4 9.34 823 351 0.0274 38.0 3.8 34.2 
102 82.6 11.3 816 342 551 0.0332 37.4 Onl 33.7 
103 104.6 8.97 858 340 596 0.0263 34.0 4.0 30.0 
105 30.4 30.8 846 302 609 0.0903 29.2 3.7 25.5 
106 40.0 23.4 $39 302 602 0.0685 29.6 Sul 25.9 
107 28.8 32. 835 300 ° 598 0.0955 29.5 Sd 25.8 
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TasBLe 7 (CoNnTINUED) 
Heat TRANSFER TO CLoups or SAND PARTICLES 


o Rm | hy | hm — hy 
6 . tw A At NA 
Run lb./(min.) Im N Ap 
sec. deg. F. | deg. F. | deg. F. pedis 
(sq. ft.) | es g Cee ae aye B.t.u./(hr.) (sq. ft. particle 
surface) (deg. F.) 
b. Particle size—40-60 mesh; 100 gm. fed 
108 32.6 28.8 839 339 581 0.1037 28.4 3.8 24.6 
108b 36.4 25.8 832 340 572 0.0929 28.8 Soene 25.1 
109 46.3 20.2 827 360 554 0.0727 32.0 3.8 28.2 
110 121.4 Ut 824 414 514 0.0278 40.6 3.7 36.9 
111 101.9 9.21 868 239 545 0.0332 40.4 4.4 36.0 
112 191.0 4.91 860 -| 468 517 0.0177 46.1 4.5 41.6 
114 102.4 9.15 848 417 540 0.0330 38.8 4.2 34.6 
JS) 85.5 10.9 $43 422 548 0.0394 38.6 4.1 34.5 
116 53.4 17.6 859 378 572 0.0634 32.8 4.1 280 
c. Particle size—60—-80 mesh; 100 gm. fed 
ah We ane 12.9 853 451 spall 0.0961 21.2 4.2 a0 
118 Od 12.2 842 452 518 0.0924 21.1 4.2 9 
118b 173.6 5.40 841 505 486 0.0402 26.8 4.5 3 
119 200.4 4.67 841 532 447 0.0348 29.9 4.7 2 
120 150.5 6.23 8538 540 455 0.0464 29.9 4.8 mL 
121 116.4 8.05 848 501 481 0.0600 26.0 4.5 a) 
122 43.8 21.4 842 428 526 0.1596 19.8 4.1 He 
123 114.0 8.23 833 499 466 0.0614 26.6 4.5 sal 
124 61.6 15.2 824 439 498 ORs 21.4 4.0 4 
125 Poe) 33.6 826 351 558 0.250 14.9 B76 2 
C. Small Furnace; Nominal Wall Temperature 700 deg. F. 
a. Particle size—30—40 mesh; 100 gm. fed 
180 359 .2 2.61 (2a 346 455 0.0076 45.7 3.2 42.5 
181 142.0 6.60 718 318 467 0.0193 40.4 Saul 37.4 
182 250.0 Ouho: 716 332 456 0.0110 43.7 oil 40.6 
183 99.6 9.42 714 300 475 0.0276 37.3 Bydl 34.2 
184 62.2 1a 715 286 485 0.0442 34.4 3.0 31.4 
185 32.0 29.3 714 249 506 0.0858 28.2 2.9 25.3 
186 28.6 32.8 714 246 507 0.0960 27.9 Pets} 25.1 
195 51.0 18.4 698 200 484 0.05389 30.6 2.8 27.8 
b. Particle size—40-60 mesh; 100 gm. fed 
| 
170 58.7 16.0 698 318 447 0.0575 35.0 3.0 32.0 
bgt 89.8 10.5 695 340 426 0.0378 38.6 3.1 35. 5 
172 59.7 lta, 7 6938 310 446 0.0565 33.4 2.9 30.5 
173 66.8 14.0 693 320 438 0.0504 BOLO 3.0 32.6 
174 150.0 6225: 695 384 412 0.0225 43.5 Bac) 40.2 
175 ole. 3.0 693 382 396 0.0108 48.0 Ome 44.8 
L7G 216.7 4.35 691 SMS) 399 0.0157 46.5 3.2 43, * 
iW ety 112.8 8.32 697 342 429 0.0299 aha il 3.0 36.0 
178 30.5 30.7 697 260 482 0.1105 25.0 | 2.8 22. 7 
179 33.8 Pi ens) 698 280 270 0.1000 ZSidt |) 229 25.4 
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TaBLE 7 (CONTINUED) « 
Heat TRANSFER TO CLOUDS oF SAND PARTICLES 


Pe e Ieee | hy | hm eae hy 
6 : é A Alim D 

Run lb./(min.) i 

sec. ee ft.) deg. F. | deg. F. | deg. F. sq. ft. Bou Greate parcels 

surface) (deg. F.) 
— rw 
c. Particle size—60-80 mesh; 100 gm. fed 
187 224.0 4.21 743 460 395 0.0314 28.6 3.8 24.8 
188 89.6 10.5 137 418 420 0.0783 23.9 3.4 20.4 
189 222.0 4.22 732 451 389 0.0314 28.5 3.8 24.7 
190 39.7 23.6 728 347 457 0.176 17.9 3.1 14.8 
191 36.4 25.8 724 320 472 0.192 Pos 3.1 12.7 
192 26.8 35.0 709 324 454 0.261 16.6 229 asi 
193 36.0 26.0 706 336 442 0.194 18.0 3.0 15.0 
194 70.0 13.4 702 368 417 0.0998 21.0 3.2 17.8 
D. Small Tube Furnace; Nominal Wall Temperature 575 deg. F. 
Particle size—40-60 mesh; 100 gm. fed 

155 28.0 83m: 586 218 392 0.1207 2508 Dae 23.5 
156 92.4 TORT 584 298 340 0.0366 42.0 2.3 39.7 
157 50.4 18.6 586 260 366 0.0670 33.6 2.2 31.4 
159 202.8 4.62 578 308 326 0.0166 45.5 2.3 43.2 
160 be? 221 578 278 346 0.0436 38.2 2.2 36.0 
161 38.9 24.1 576 239 368 0.0868 30.3 2.2 28.1 
162 12.3 13.0 579 274 350 0.0468 37.0 2.2 34.8 
163 S125 29.8 582 215 390 0.1073 20.2 Perel} 23.1) 
164 70.9 132 576 270 350 0.0475 36.6 2.2 34.4 
165 36.1 26.0 570 220 376 0.0936 Pee Pra i 25.1 
166 55.8 16.8 566 251 352 0.0605 33.4 2a 31.3 
167 128.8 7.29 563 292 320 0.0262 43.6 2/2: 41.4 
168 89.8 10.4 563 272 334 0.0374 38.8 2.2 36.6 
169 176.2 Doe 561 294 317 0.0192 44.3 2.2 42.1 


ratio, 2k/hepD,, vary from 94 to 246. Therefore, the assumption of 
negligible temperature gradient through the particle is satisfactory for 
the particle sizes studied in this investigation. 

The validity of the assumption that the overall rate of heat transfer 
is proportional to the logarithmic mean temperature difference between 
the wall and the particle for a particle falling with accelerated motion 
through a tube at constant wall temperature may be proven as follows: 

The heat balance on a single particle requires 


wpcdt = hnA» (tw — ty) dé. (63) 


If t, is constant, this equation may be integrated between limits to 
give 
5 


bes 7 tpi hwA» 
In = A6;. (64) 


tw —< typo W pC 
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TaBLE 7 (CONTINUED) 
Heat TRANSFER TO CLoups oF SAND PARTICLES 


| 
Rin | h hn — h 
° F Ie m T 
Run ree! u as a an A on ae Ne 
sec. min. eg. F. | deg. F. | deg. F. . ft. 
ieee (sq. ft.) . @ . HS B.t.u./(hr.) (sq. ft. particle 


surface) (deg. F.) 


E. Large Furnace; Nominal Wall Temperature 1000 deg. F. 
a. Particle size—30-40 mesh; varying amounts fed as indicated 


13 448 30.0 25.8 992 492 614 O.712 29.4 Dao) 23.9 
14 378 45.0 14.5 992 559 567 0.400 36.8 5.8 31.0 
15 434 90.0 8.32 981 600 524 0.230 43.2 (aul BYfaul 
16 334 120.0 4.80 973 632 492 OL1s27 49.0 Gime 42.9 
ile 443 25.0 30.5 =| 984 467 622 0.842 27.4 5.4 22.0 
18 342 180.0 3.28 997 672 486 0.445 37.1 5.9 31.2 
b. Particle size—40-60 mesh; varying amounts fed as indicated 
20 452 25.0 ol..2 985 508 596 0.1155 23.4 5.5 17.9 
2. 406 45.0 15.6 995 623 522 OL577 33.8 6.3 2020 
22 328 180.0 3.14 998 745 420 0.1161 52.3 Tath 44.6 
23 356 25.0 24.6 983 547 566 0.910 26.8 ewe a 
24 401 50.0 13.8 967 625 490 (0) tapi 36.2 6.1 30.1 
25 278 60.0 8.00 977 700. 437 0.296 46.4 6.9 39.5 
26 233 90.0 4.46 984 746 402 0.1652 54.2 stl 46.5 
F. Small Furnace; Nominal Wall Temperature 1050 deg. F. 
a. Particle size—30-40 mesh; varying amounts fed as indicated 
196 18.8 30.0 5.88 1050 484 706 0.0168 44.4 LiyaYA 38.8 
197 2875 60.0 4.45 1045 533 670 0.0127 62.5 5.8 46.7 
198 (hoses) 120.0 5.90 1044 508 685 0.0169 48.1 Sail 42.4 
199 102.8 180.0 i eeho) 1041 513 680 0.0153 49.2 ayers 43.5 
200 180.7 240.0 7.05 1039 489 692 0.0202 45.7 5.6 40.1 
201 41.2 30.0 12.9 1045 452 722 0.0369 40.0 5.6 34.4 
202 81.6 60.0 12.8 1043 456 720 0.0366 40.4 5.6 34.9 
203 116.0 90.0 12.1 1042 460 714 0.0346 41.3 5.5 35.8 
204 165.9 120.0 13.0 1037 448 715 0.0372 39.8 5.4 34.4 
205 242.7 180.0 12.7 1031 448 719 0.0363 39.7 5.5 34.2 
b. Particle size—40-60 mesh; varying amounts fed as indicated 
207 15.1 15 9.44 1035 553 644 0.0330 47.6 5.8 41.8 
208 31.2 30 9.76 1035 550 647 0.0341. 46.0 5.8 40.2 
209 43.9 45 9.15 1034 562 6388 0.0320 47.8 5.8 42.0 
210 68.5 75 8.56 1031 548 644 0.0300 45.7 Bnd 40.0 
| 
G. Small Furnace; Nominal Wall Temperature 725 deg. F. 
Particle size—40-60 mesh; varying amounts fed as indicated 
| 
228 255 60.0 39.9 718 268 502 0.140 26.1 2S 
229 152 40.0 35.6 Chie 279 490 0.125 2E8 220 
230 84.5 20:2 39.2 719. Pal 500 0.137 26.6 2.9 
231 42.7 10.0 40.1 719 274 495 0.140 27.0 2.9 
232 48.9 120.0 3.82 (OR) 395 418 0.0134 48.1 3.2 
Zoo. 90.5 240.0 OEOe 721 408 409 0.0124 51.3 Suz 
234 8.9 30.0 2.78 722 411 407 0.0097 52.1 See 
235 37.4 100.0 3.51 721 409 408 0.0123 51.6 3.2 
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TaBLE 7 (CONCLUDED) % 
Herat TRANSFER TO CLOUDS OF SAND PARTICLES 


FP hn | hy | hm — hr 
Mca ee b./ bu A_ | Atm | NA» 
“a. sec. ae deg. F. | deg. F.| deg. F. | sq. ft. Pee CONC eS 


surface) (deg. F 


H. Small Furnace; Nominal Wall Temperature 1050 deg. F. 
Particle size—40-60 mesh; varying amount fed as indicated; atmosphere of CO2 


211 89.9 75 ibe 1035 576 629 0.0427 45.9 622 39.7 
212 35.9 45 7.48 1038 592 621 0.0282 48.2 6nd ehh 
213 9.5 15 5.94 1041 649 582 0.0224 57.5 On7 50.8 
214 Palas 30 6.65 1042 637 593 0.0251 54.9 6.7 48.2 
215 94.4 120 7.38 1034 617 599 0.0279 523 (Sof 45.8 
1 
I. Small Furnace; Nominal Wall Temperature 725 deg. F. 

Particle size—40-60 mesh; varying amounts fed as indicated; atmosphere of CO2 
39.6 10 37.8 732 292 497 0.1407 26.8 Soul 23).7 
ficval 20 36.3 730 302 490 0.1350 26.8 2.9 23.9 
164.4 40 39.2 727 288 495 0.1460 26.4 3.0 23.4 
240.3 60 38.2 723 283 495 0.1421 26.0 3.0 23.0 


Here Ad; is the time of fall. The true average temperature difference 
is seen to be the logarithmic mean if the expression for the total heat 
transferred in time is written 


LO (Eo) ba) eA eth NO (65) 


By comparison with Equation (64) it is seen that 


ayy tpl 
(S a bo) ae = sori ear a (as a aS veo (66) 
es 
Ubon tp 


The derivation of Equation (64) is based on the assumption of a 
constant value of hm». Numerical integration of Equation (63) using 
for h» theoretical values obtained from individual coefficients calcu- 
_ lated by means of Equations (36), (57), and (58) showed that under 
“ty pical conditions h», might increase from 30 to 36 B.t.u./(hr.) (sq. ft.) 

(deg. F.), and remain constant at the latter value over the lower 
half of the furnace. This change was too small to affect the validity of 
the logarithmic mean assumption. 
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TABLE 8 
Hear TRANSFER TO CLOUDS or CARBORUNDUM PARTICLES 


Nn | hy | Rm — hy 
Run Oe eae Hote e iaee es eee 
sec. : ; eg. F. | deg. F. | deg. F. sq. ft. 
(sq. ft.) : bs . a B.t.u./(hr.) (sq. ft. particle 
surface) (deg. F.) 
A. Small Furnace; Nominal Wall Temperature 1050 deg. F. 
a. Particle size—No. 40; 100 gm. fed 
C19 496.0 1.89 1052 627 611 0.0069 55.1 6.8 48.3 
C20 193.0 4.86 1046 577 641 0.0179 47.4 6.4 41.0 
C21 116.2 8.09 1039 539 658 0.0296 42.9 6.2 36.7 
C22 72.2 13.0 1038 506 674 0.0475 39.0 5.9 33.1 
C23 58.2 16.1 1032 498 676 0.0593 38.0 5.8 32.2 
C24 38.0 24.7 1020 452 692 0.0905 33.2 ye) PHP 
25 49.2 19.0 1017 473 675 0.0700 35.8 5.6 30.2 
C26 29.0 32.4 1020 429 706 0.1185 30.6 5.4 25.2 
b. Particle size—No. 60; 100 gm. fed 
C27 286.0 3.28 1047 697 bos 0.02388 3.4 7.2 26.2 
C28 146.4 6.40 1035 662 565 0.0465 2:0 6.9 Zoe 
C29 118.0 7.95 1028 641 578 0.0576 2 6.6 23.6 
C30 78.0 12.0 1022 613 590 0.0870 AG ‘Ono 22.0 
C31 59.2 15.8 1021 583 586 0.1146 315) 6.1 19.4 
C32 42.0 22).3: 1022 548 635 0.1618 sth 6.0 16.7 
C33 384.9 28.9 1026 543 643 0.1950 2.4 5.9 Gms: 
C34 27.9 33.6 1028 512 666 0.244 2 5.7 14.5 
B. Small Tube Furnace; Nominal Wall Temperature 725 dee. F. 
a. Particle size—No. 40; 100 gm. fed 
Cl 45.2 20.8 720 290 475 3.8 310) 23.8 
G2 21.8 43.0 714 283 495 ao, 2.8 20.7 
C3 27.3 24.4 707 289 486 F .8 2.9 21.9 
C4 39). 7 23.6 705 282 475 0.0869 5.6 2.9 23.7 
(O13 273.0 3.44 709 375 420 0.0127 m0) 3.3 SOnd 
C6 91.0 15.4 705 317 454 0.0566 9 3.0 28.9 
C7 93.7 10.0 745 364 464 0.0368 520) Sey ol.7 
C8 180.6 6.19 743 393 443 0.0191 me) on) 38.4 
C9 155.0 6.05 738 385 442 0.0222 5 3.5 38.0 
b. Particle size—No. 60; 100 gm. fed 
C10 207.8 4.52 (ey 463 380 0.0328 29),2 3.0 25.95 
Ci OUOn a: 2.50 728 471 368 0.0181 30.4 Sl al een 
C12 124.8 Venray) 125, 434 894 0.0545 26.4 B.D 22.9 
C13 97.3 9.65 G2 415 403 0.0700 24.5 3.4 21.1 
Cl4 70.2 13.4 712 385 414 0.0970 PPX A Bee 18.9 
C15 46.2 20.3 a? 358 434 0.1472 18.5 Bil 15.4 
C16 26.4 33.0 710 320 455 0.239 16.4 2.9 13.5 
C17 20.8 45.1 707 315 488 0.327 ilsjqal 228 | 12.3 
C18 39.6 24.3 707 343 440 0.1762 18.4 (330) 15.4 


12. Analysis of Data—When a cloud of particles falls through < 
heated gas space under the influence of gravity, the time of contact of 
a single particle with the hot gas is essentially independent of the rate 
of feed. On the other hand, the ratio of the area of the particles to 
the area of the wall, y, should be directly proportional to the feed rate. 
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TABLE 9 
Heat TRANSFER TO CLOUDS or ALOXITE PARTICLES 


‘ 


hin | hy | hm ae hy 
Run @- lip./(min] 98 Deon) Meee 
sec. Ga: ft.) deg. F. | deg. F. | deg. F. sq. ft. Bean) (eae pareicls 
surface) (deg. F.) 
A. Small Furnace; Nominal Wall Temperature 1050 deg. F. 
a. Particle size—No. 46; 100 gm. fed 
Al6 395.0 2.38 1032 563 634 0.0076 52.9 (oy. J1 46.4 
Al7 187.2 5.62 1031 525 660 0.0179 47.5 5.9 41.6 
A18 92.6 10.1 1027 454 702 0.0323 37.8 5.6 31.2 
A19 214.8 4.36 1073 550 692 0.0140 50.5 6.5 44.0 
A20 151.0 6.21 1072 537 693 0.0199 48.3 6.4 41.9 
A21 103.8 9.05 1070 508 713 0.0290 44.1 6.2 37.9 
A22 68.8 13.6 1070 467 740 0.0435 39.3 6.0 33.3 
A23 46.3 20).2 1072 425 769 0.0646 34.0 5.8 28.2 
A24 36.2 25.9 1070 400 780 0.0829 31.1 ert 25.4 
A25 30.2 31.0 1062 380 773 0.0992 29.4 5.5 23.9 
A26 24.2 38.8 1082 356 800 0.1240 26.4 5.4 21.0 
b. Particie size—No. 60; 100 gm. fed 
A27 197.2 4.75 1043 613 613 0.0220 44.0 6.6 37.4 
A28 252.0 Sele 1038 631 593 0.0172 48.4 6.9 41.5 
A29 90.6 10.3 1069 555 681 0.0524 32.6 6.4 26.2 
A30 54.2 ef 363 1062 508 704 0.0799 31.2 Gal Zoe 
A31 64.0 14.6 1055 527 685 0.0679 by) 6.1 Pia 
A32 35.6 26.3 1051 453 726 0.1214 31.8 Gal 25.7 
A33 31.6 29.7 1043 425 737 0.1371 29.3 5.8 23.5 
B. Small Furnace; Nominal Wall Temperature 700 deg. F. 
a. Particle size—No. 46; 100 gm. fed 
A8& 289.2 3.24 716 340 450 0.0104 43.5 3.2 40.3 
AQ 111.6 8.41 714 306 470 0.0269 BY oe) 3.1 34.1 
Al0 184.0 5.10 713 320 461 0.0163 40.6 3.2 37.4 
All 67.4 13.9 az 280 483 0.0445 32.8 eno) 29.8 
A12 32.8 28.6 711 230 513 0.0915 27.6 2.8 24.8 
A13 43.6 21.4 711 247 502 0.0688 27.8 29 24.9 
Al4 21.5 43.6 712 205 527 0.1396 21.4 2.7 18.7 
Al5 17.6 5320 715 196 530 0.1704 20.5 256 17.8 
b. Particle size—No. 60; 100 gm. fed 
Al 178.0 5.26 711 382 418 0.0244 36.6 3.3 33.0 
A2 225.6 4.15 704 393 403 0.0192 39.4 3.3 36.1 
A3 74.0 120 698 334 435 0.0582 30.4 3.0 27.4 
A4 103.2 9.08 692 354 417 0.0420 33.4 ‘On 30.3 
Ad 29.2 ol 687 270 465 0.1483 22.8 2.8 20.0 
A6 36.8 PASE) CP 288 490 0.1179 23.1 3.0 20.1 
A7 51.0 18.8 718 313 470 0.0869 25.2 3.0 22.2 


Now if we may assume as an approximation that the individual coef- 


ficients, Nep and he», are essentially independent 


of the feed rate, it 


follows that the ratio of the convective resistance at the wall to that 
at the particle, NA,hcp/Awhew, must likewise be proportional to the 
feed rate. As shown in Tables 7, 8, and 9 for each series of runs the 
feed rate was varied over a range of nearly fifteen-fold while other 
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conditions were maintained constant. The corresponding variation in 
the resistance ratio thus permits estimation of the individual convec- 
tion coefficients. For this purpose either of the following equations 
might be used: 


ae 67 
Ve =x Ny i: hs Now \ 

1 1 1 
= 2 (68) 


oy Wes. — h,) Myer Vhep , 


For the first equation a plot of 1/(h» —h,) against F is required, 
from which the reciprocals of h-, and hex, respectively, may be deter- 
mined from the slope and the intercept at F = 0. For the second equa- 
tion, 1/y(hm — h,) is plotted against 1/F. The slope of the line now 
gives 1/h-», and extrapolation to the zero abscissa, representing infinite 
feed rate, gives 1/Acw. It should be noted that in this equation 
y(hm» —h,) represents the overall convection coefficient based on the 
wall area. Consequently, any error in computing the time of fall and 
area of the particles will not affect the value of he» obtained by the 
extrapolation. Possibly for this reason, and also because the assump- 
tion of the independence of the coefficients of feed rate apparently holds 
better at high feed rates, the second method of plotting the data was 
more satisfactory than the first. Before examining the data it is well to 
consider the validity of the foregoing assumption as far as possible 
from the theoretical standpoint. 

According to Equation (36), which is valid for the range of 
Reynolds numbers covered in the experimental work, he» is affected by 
the feed rate indirectly through the temperature of the gas, since for 


48 ILLINOIS ENGINEERING EXPERIMENT STATION 


Small Furrrace — Large Furnace- 

A -4,=/050°F., Air V~ ty = 1000, Air 
My ODO COs 
D-#4,= 850°F, Al/r 
ot, = 700°F, Air 
e -4,= 700°F, Coz 
= n= UO day Ce 


Wag lad) 


60-80 Mesh 


Fic. 9. GRAPHICAL REPRESENTATION OF Herat TRANSFER 
DaTA FOR SAND PARTICLES 


any wall temperature the gas temperature is decreased as the feed rate 
is increased. The effect, however, is negligible since, for air, the group 
V kVpc/D, varies only about one per cent for each one hundred 
degrees Fahrenheit. 

In the case of hey the effect of the feed rate is not so obvious. If 
only natural convection from the wall is considered, for a constant wall 
, temperature, the Grashof group is expected to increase as the feed rate 
is increased. The approximate extent of this effect is shown in Fig. 8 
in which the calculated values of the coefficient are plotted against the 
feed rate. For the conditions chosen the relationship is seen to be 
nearly exponential. The slope of the logarithmic line, however, is only 
0.15, so that, from this standpoint, the approximation is nearly true. 

The graphical representation of the complete data according to 


HEAT TRANSFER TO CLOUDS OF FALLING PARTICLES 49 


700°F., Small Furnace 
%, Small Furnace 


Fig. 10. GrapHicAL REPRESENTATION OF Herat TRANSFER 
Data FoR CARBORUNDUM ParTICLES 


Equation (68) is shown in Figs. 9, 10, and 11. It will be observed 
that, for a given furnace, the points fall very near a single straight 
line regardless of the wall temperature. The change in slope from one 
furnace to another is due principally to the change in the ratio between 
y and F. Using the method of Least Squares for the linear equation 
the coefficients, he» and hep, were calculated for each of the series of 
runs. The results are shown in Table 10. It is apparent that neither of 
the coefficients varies with the wall temperature beyond the experi- 
mental error. Consequently, average values are justified for each 
particle size. 

In Table 10 the values of hey and h,, calculated from Equations 
(58) and (35), respectively, are also given. The observed values of Rew 
range from 2.1 to 3.9 B.t.u./(hr.) (sq. ft) (deg. F.) for the small 
furnace and from 3.8 to 4.1 for the large furnace. The average coeffi- 
cient of natural convection calculated from the Rice equation is 
about 1.0, and is independent of the diameter of the furnace. Undoubt- 
edly this discrepancy is due to the increased convection caused by the 
particles themselves, and by the cooling of the gas at the top of the 
furnace. In a closed furnace there must be a circulation of the gas up the 
walls and down the center of the furnace. The convection at the walls 
apparently increases with the size of the furnace. Unfortunately, two 
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furnace sizes are not sufficient to establish even an empirical relation- 
ship for this effect. For values of the Grashof group between 10* and 
10°, a rough approximation might be made by multiplying the coeffi- 
cient for natural convection by a factor equal to the diameter of the 
furnace in inches. For values of the group greater than 10° it is 
probably better to omit the correction factor entirely, and the overall 
coefficient would then be conservative. The latter range includes most 
cases of practical importance for which the furnace is large and heat 
transfer by radiation predominates, so that even large errors in the 
calculated value of he» can be neglected. Furthermore, in many 
furnaces heat is not transmitted through the walls but is introduced 
directly in hot gases so that he» does not enter into the calculations. 
The observed values of hep agree within 20 per cent with those 
calculated from Equation (35) for all cases except for the smallest 
particles where they are about 50 per cent low. Since the observed 
values depend directly on the calculated time of fall, a possible expla- 
nation of the discrepancy lies in the downward motion of the gases at 
the center of the tube which would result in a shorter time of contact 
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TABLE 10 
SuMMARY or Resutts on Inpivipuan Hear Transrer CoprricreNTs 
Rew, Bete: hep, B.t.U. 
’ Wall (hr.) (sq. ft. wall area) (hr.) (sg. ft. particle 
Material pmern Temperature (deg. F.) area) (deg. F.) 
ize 
deg. F. 
Calculated Calculated 
Observed (Equation Observed (Equation 
(58)) (35)) 
A. Small Furnace 
1050 2.4 0.64 to 0.98 48 4 
Sand 30-40 850 2.3 0.70 to 1.02 48 33 
700 2.3 0.63 to 1.01 47 53 
Average 2.3 48 
1050 2.5 0.75 to 1.04 52 59 
Sand 40-60 850 2.8 0.74 to 1.02 49 59 
700 2.8 0.69 to 1.05 55 58 
575 3.4 0.75 to 1.04 52 58 
Average 2.8 52 
: 1050 2.9 0.72) to: 1.10 30 67 
Sand 60-80 850. 2.6 0.76) to 1.05 34 66 
700 3.0 73 to. 1,05 32 65 
Average 2.8 32 
1050 3.5 0.61 to 1.02 44 55 
Carborundum No. 40 700 3.0 0.72 to 1.08 46 5 
Average 3.2 45 
1050 3.9 0.69 to 1.02 30 65 
Carborundum No. 60 700 $0) 0.65 to 1.04 27 63 
Average 3.4 | 28 
: 1050 Papal 0.63 to 1.05 53 57 
Aloxite No. 46 700 2.5 0.68 to 1.10 45 56 
Average DO) 49 
1050 2.4 0.72 to 1.03 47 64 
Aloxite No. 60 700 ne] 0.72 to 1.06 51 62 
Average 238 49 a 
B. Large Furnace 
Sand 30-40 1000 4.1 O73 to: L103 49 54 
Sand 40-60 1000 3.8 0.75 to 1.04 52 54 


than that calculated. The error would be greater for the smaller 
particles than for the larger ones. It is unlikely that an error in the 
measurement of the particle areas would be positive, thus resulting 
in small coefficients, since concave surfaces are known to exist in all 
of the particles to some extent. 

The values of h,, for carbon dioxide at a wall temperature of 
1050 deg. F. are definitely greater than those obtained for air, as shown 
in Fig. 9, but the values at 700 deg. F. are essentially the same for 
both gases. If it is assumed that the radiation interchange between the 
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eas and the suspended particles is negligible,,and therefore that hep is 
the same for the two gases, the coefficient of heat transfer for radia- 
tion from the furnace wall to the carbon dioxide can be estimated by 
means of the usual equations for gas radiation. The values obtained 
in this way for the two temperatures are, respectively, 1.0 and 0.25 
B.t.u./(hr.) (sq. ft. wall area) (deg. F.). The higher value would 
represent an increase of about 40 per cent in the value of he» obtained 
without superimposed radiation. While the data are insufficient to 
determine the exact values of the individual coefficients by plotting, 
it appears that the larger values of h,, obtained when carbon dioxide 
is the ambient gas may be explained on the basis of the radiant inter- 
change between the wall and the gas. 


V. CoNcLUSIONS 


13. Summarized Conclusions—The foregoing discussion may be 
summarized as follows: 

(1) The rate of heat transmission from a hot furnace wall to 
clouds of falling particles is rapid. When the diameter of the furnace 
is small, the transmission takes place principally by convection from 
the wall to the gas and from the gas to the particles. The fraction 
transmitted by radiation increases rapidly with the size of the furnace 
and the temperature of the wall (q. infra). The factor producing the 
most important effect on the overall transmission is the rate at which 
the particles are fed to the furnace. When radiation is small, the main 
effect of varying the feed rate is to change the ratio of particle area 
to wall area, thus changing the relative magnitudes of the thermal 
resistances at the wall and at the particle without altering appreciably 
the respective convection coefficients themselves. 

(2) The coefficient of heat transfer by convection from the gas 
to the particles is a function of the particle size and the rate of fall, 
and is nearly independent of the temperature of the ambient gas. This 
coefficient may be calculated from fundamental considerations of heat 
transfer to a submerged body in a moving fluid. The general equation is 


k co [= enn,” 
lise = SS 
D; i=0 1 + b4 


(35) 


where the constants b; are defined by the relationship: 2b;2Re:Pr = 
the squares of the roots of the first order Bessel function. For particles 
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having a specific gravity of about 3, and larger than 400 microns, fall- 
ing through air, the foregoing equation may be simplified to 


| kV pc 
hep = 0.714 a (36) 


Pp 


(3) The coefficient of heat transfer by convection from the furnace 
wall to the gas is several times larger than the values predicted for 
natural convection alone, and is apparently a function of the size of 
the furnace. For furnaces less than eight inches in diameter and at 
wall temperatures less than 1200 deg. F., it is recommended that the 
natural convection coefficient be multiplied by a constant equal to the 
diameter of the furnace in inches. For larger furnaces, or for higher 
wall temperatures, the heat transmission takes place primarily by 
radiation, and the correction may be neglected. 

(4) The radiation coefficient is a function of the concentration of 
particles in the cloud and of the shape of the furnace. The coefficient 
may be calculated from the equation 


0.172e, [(7.,/100)4 — (T,/100)*] 


57 
af Ch: ve ih) 


r 


Values of the absorptivity of the cloud, «,, are given in Table 1 for 
spherical, cylindrical, and rectangular furnaces. 


14. Design of Large Furnace ——In order to illustrate the applica- 
tion of the foregoing equations an example will be taken from the 
previous publication.* The problem is the design of a flash calciner 
for decomposing 83.2 tons of zinc sulphite hydrate per day into zinc 
oxide, water vapor, and sulphur dioxide. The conditions are as follows: 

Feed, dry ZnSO;:21%4H.O enters at 70 deg. F. 

Particle size, 200 microns, found by experiment to be the arith- 

metic mean of particles just passing 100 mesh standard screen. 

Dehydration temperature, 200 deg. F. 

Decomposition temperature, 500 deg. F. 

Specific heat of feed, 0.24 B.t.u./(1b.) (deg. F.). 

Heat of dehydration, 40 800 B.t.u./1b. mole. 

Heat of decomposition, 58 100 B.t.u./1b. mole. 

Emissivity of particles, 0.5. 


*Univ. of Ill. Eng. Exp. Sta. Bul. 324. 
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Fic. 12. Hear TraNsrer Corrricients aNp CaLcuLATED Hericuts or 
RECTANGULAR FURNACES FOR CALCINATION or ZINC SULPHITB 


A feed rate of 11.7 Ibs. per min. per sq. ft. is assumed so that the 

, total cross-sectional area of the furnace will be 10 sq. ft. The furnace 
is to be constructed in the form of rectangular retorts with walls of 
high heat conductivity. Silicon carbide tile are available, 1144 inch 
thick, tongue and groove, with a heat conductivity of 105.5 B.t.u./(hr.) 
(sq. ft.) (deg. F./in.) at 1110 deg. F. and 112.5 B.t.u./(hr.) (sq. ft.) 
(deg. F./in.) at 2900 deg. F. The retorts may be placed in a row with 
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the common walls formed of slotted tile. The walls will be sealed by 
compression with suitably spaced springs. The furnace will be heated 
by horizontal gas muffles along each side, the flames traveling in several 
passes from bottom to top. The outside walls will be constructed of 
fire-brick and will be insulated. 

Calculations were made to determine the effect of the shape and 
size of the retorts on the height required for a heat input sufficient to 
completely dehydrate and decompose the particles. These calculations 
were made by applying Equations (5), (57), and (58) to successive 
short sections, taking into consideration the change in the velocity, 
mass, and temperature of the particles and of the gas. The tempera- 
ture of the furnace wall was considered to be constant. The evolved 
gases were assumed to leave at the bottom of the furnace so that there 
was an increase in the mass velocity of the gases in the two reaction 
zones. The reactions were considered to take place only when the 
particles reached the temperatures at which the respective decompo- 
sition pressures are one atmosphere, and at these temperatures the 
rate of reaction depended only on the rate of heat input. 

The results of these calculations are plotted in Fig. 12. The fol- 
lowing conclusions may be drawn from this study. 

(1) With a uniform wall temperature of 1000 deg. F. the height of 
the calciner of rectangular cross section 30 inches long and 12 to 20 
inches wide approaches 40 feet. Nothing is to be gained by using 
smaller cross sections, because the length does not decrease as rapidly 
as the area of the cross section decreases and, therefore, the increased 
number of units required more than offsets the decrease in height. 

(2) An increase in wall temperature to 1200 deg. F. decreases the 
height by nearly 50 per cent. 

(3) In these large rectangular furnaces the heat transfer is ninety 
per cent by radiation. If the source of heat is located at the bottom, 
a ereater heat flux can be obtained than is shown by the curves, 
because of the great effect of the temperature level on radiant heat 
transfer. Furthermore, since the greatest part of the heat is required 
at the bottom in the decomposition zone, upward flow of the combus- 
tion gases is highly desirable. 
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